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ABSTRACT 

The purpose of this Investigation is to evaluate imaging radar as 
an adjunct to conventional petroleum exploration techniques, especially 
linear mapping. Linear features were mapped from several remote sensor 
data sources including stereo photography, enhanced Landsat imagery, SLAR 
radar imagery, enhanced SAR radar Imagery, and SAR radar/Landsat combin- 
ations. Linear feature maps were compared with surface joint data, 
subsurface and geophysical data, and gas production in the Arkansas part 
of the Arkdiiia basin. 

The best Landsat enhanced product for linear detection is found 
to be a winter scene, band 7, uniform distribution stretch. . Of the 
•individual sAR data products, the VH (cross-polarized) SAR radar mosaic 
provides for detection of most linears*, however, none of the SAR enhance- 
ments is significantly better than the others, Radar/Landsat merges 
may provide better linear detection than a single-sensor mapping mode, 
but because of operator variability the results are inconclusive. 
Radar/Landsat combinations appear promising as an optimum linear mapping 
technique, if the advantages and disadvantages of each remote sensor are 
considered. 

Although interpreters generally agreed on the or.ientation of major 
linear trends i they differed greatly on other parameters, such as total 
number, average length, and coincidence of linears. Mapping of relative 
concentrations of linears (density) appears to be somewhat less influenced 
by operator variability. 

In the southern region of the Arkansas Arkoma basin, little 
relationship is discernible between surface structure and gas production, 
and no correlation is found betv;een gas production and linear proximity 
or linear density as determined from remote sensor data. However, in 
the northern Arkoma basin, a positive correlation is found between linear 
density and gas production from a high-chert-content carbonate reservoir. 
Linear density analysis should not be the single basis for petroleum 
exploration on the northern flank of the Arkoma basin. Regional and local 
structural position and surface joint orientation are equally important 
data sets to be evaluated in an exploration program. However, the proba- 
bility of economic savings in the exploration for fractured reservoirs 
in other parts of the work can be significant if linear density maps, 
derived from radar and Landsat merges, are used to target prospects for 
furthor geological and geophysical investigations. 
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I. INfROI^lJCTION 

The space program has oonflcmed that most remote sensing 
techniques previously proved useful from aircraft altitudes 
have equal or greater value from orbital altitudes, Landsat 
imaging systems have been used to demonstrate that earth 
observations from space can provide valuable information for 
the exploration and management of energy and minefal resources. 
However# even a brief review of the present state-of-the-art 
i of earth resource exploration from space shows that an 
improvement must be made in spaoeborne observing and measuring 

. f: 

techniques. Present remote sensor technology for geologic 
exploration is simply inadequate to aid in the discovery of 
sufficient new reserves to satisfy the accelerating demand 
for petroleum and natural gas. NASA has recognized this 
fact and several satellites and space mission’s are being 
readied for launch during the next few years to take advantage 

j. 

of new sensor technology. New sensors are being designed 
that are adapted more specifically to the earth resources 
applications. Spaoeborne microwave systems appear to have 
immediate application for geologic mapping. 

During the past two decades the petroleum industry has 
shown increasing interest in fractured reservoirs because 
many of the world's oil and gas fields have reservoirs in 
which porosity and# perhaps more important, permeability 
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are fracture controlled. Knowledge of fracture control in 
a reservoir not only aids the primary recovery of hydrocarbons 
but also guides the design of well stimulation and enhanced 
recovery programs . Because fractures are commonly 
propagated upward and reflected at the earth's surface as 
subtle linearSf detection of these surface features is 
extremely important in many phases of petroleum exploration 
and development. 

As an adjunct to field investigatibns, various remote 
sensor technigues have been used for fracture analysis and 
geologic reconnaissance surveys. Photography has been the 
most useful survey system because its application has been 
refined, and a great nuiTijer of people are experienced in 
analyzing the images obtained. From the vantage point of 
the aerial camera/ geologists have been able to obtain fracture 
information that is not readily discerned from the ground. 
However/ within the last decade/ a variety of more 
sophisticated remote sensors have been used in geologic 
studies. More recently^ Landsat satellites have provided 
the impetus for renewed interest in linear tectonics. 

Very little researeh has been devoted to assessing the 
geologic linear mapping capability of radar imagers, though 
such mapping has received considerable attention in relation 
to visible and near-visible spectrum sensors . The linear 
mapping potential of multifrequency imaging radars from 
satellites appears to be very promising for petroleum 


exploration# especially in the United States. Active remote 
sensors such as multifrequency imaging radars now are heing 
recogniised as the next logical and porhaips the most 
sophisticated sensor package to be placed in space. These 
microwave remote sensors^ could have a significant impact 
on petroleum exploration in the next decade. 

To document the utility of microwave analysis for 
petroleum exploration# the Arkansas port of the Arkowa 
basin was selected as a prime test site. The research plan 
included comparison of aircraft microwave imagery , stereo 
photography# and Landsat imagery in an area where significant 
' subsurface I'orehole gecphysical data were available. The 
test site provides ocutras ting rock types and structures in 
an area where surface cover includes cropland and deciduous# 
coniferous# and mixed forest types. There is little 

discernible relationship between structure and gas production 

,•/ • 

in the area# and, stratigraphic entrapment appears to be the 

. jj ^ ^ ' 

primary factor influencing the presence of gas. 

Space imaging radar systems will provide the longest 

■ , V : : ■ _ ■ ■ ■ ■ ' ^ ^ ^ ■ ■ ■ 

wavelengths of the family of instruments used for sensing 
the electromagnetic radiation reflected or emitted from the 
earth's surface. These microwave systems should yield 
unique terrain data to facilitate geologic interpretation 
and thus will provide an important new technique for petroleum 
exploration . In anticipation of Space Shuttle programs for 
the 1980s# radar systems should be evaluated as an adjunct 
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to conventional, petroleum exploration technicjueSi especially 
linear-lineament mapping. Such an evaluation is the overall 
objective of the program described in this report. 

SubobjectivcB include: 

1) Compilabion of linear featni^es maps from appropriato 

0 

remote sensor data sets. 

2) comparison of linear features maps with surface 
joint data, subsurface structure and geophysicfal data, and 
gas productivity. 

3) Determination of the relationship between remote 
sensor linear proximity and density to gas productivity. 

4) Definition of optimum remote sensor data set for 
linear enhancement. 

5) Evaluation of the role of operator variability in 
linear feature mapping . 

II. LOCATION AND GEOLOGY OF ARKOMA BASIN 

The test site is in northwestern Arkansas, generally 
within a rectangle bounded by latitudes 35*45 'N and 35*07 'N 
and longitudes 94*30 'W and 92*45 'W, encompassing approximately 
11,000 km . It includes all of the major gas fields in the 
Arkansas part of the Arkoma basin. The Arkoma basin is an 
arcuate structural trough 400 km long extending east-west 
across central Arkansas and northeast-southwest in eastern 
Oklahoma (Fig. 1) . The basin is typified by imbricate thrust 
faults in the south and normal faults, monoclines , and 
east-west-trending anticlines and synclines in the north. 
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The Arkoma basin kanqes in width from 32 to 70 km measutod 
from notth to south. It is boundad on tho north by tho 
northern Arkansas structural platform and on the south by 

the Ouachita Mountains (Fig. 2). 

The northern Arkansas structural platform, composed of 

essentially horizontal strata of Mississippian and Pennsyl- 
vanian age, trends genorally east-west across northwestern 
Arkansas. Structural features of the platform include low 

amplitude east-west-trending folds, northeast-southwest- 
trending normal faults, and east-west-trending faults along 
the southern margin. Other prominent features of the platform 

are a series f regional northeast-southwest-trending 
topographic lineaments (aligned straight stream valley 






NO R T H E R N ARKANSAS 



S T R U C T U R A L PL AT FORM 


1000 O OO Q Q O O 0 0 0 OO OOho 000 Q C 0 0 O 0 O ^TD tJO Q*D0 OTTO OTTOTJ 0 0 00^ OIOT 

PQOooooooooodoopoopooooo >000000000000600 troooooooooopo I 
iooooqooojDOoooacfoooooooogoooooooooooooooo ^oooopoooDqoc 
Oooooooooooor/m>baixao*>^aa ooooooooooooi 

gqp 0 o 0 o cioh crot^ o o o o o o o c o o o o o o o 0*0 aao on lo o o o o o o o o olo c 

do Qi>^(TOTTa O Q do 00000 00 000 )0060600000000 ob C^O*aa^ 0 O 06 0 1 

^ m>t) 6 O O 0060 0 oHo o o o o o o o o o < o o o o o o o o o o o o o o d o 00 0 O O U 0 croorni 0 

kfco 6 O O 0060000 0 olo 0 0 O 0 0 0 O O O > 0 O o O O 0 0 o D O 0 O 0 qo 0 0 D 0 O 6 0 0 O 0 6 06 1 * 


* t* ^ w w w V w vf V 4 w VI M w V M V W wvi w w W u u V* V# v» » XV r - ' * ” 4 » y " ' v; u v w u VJ 

3 fcc>oooooocumql) 006 ooo 6 oooboooooo?)oooC arkSVI G^ooooooo 

^ 060600 00 O O j> 6.0 O^XLOP 06 6 06 0 0 0 0 6 000 0 O Q ao 0 0 6 O 0 

000060066600 600000006 OQOOOoOOOOOOOOOCQOOOOOOOQOOOOlOOt 
> C £0 O O D 0 06 O O 0 6 0 O 0 0 O O O op 6 o O O O 0 0 O 0 6 6 06006 6 00 0 0 0 O 0 O 6 6 60 0 0 
r ♦ * 0 O O op 0 o o o poo^oo o o o o o o o o o o o o o o o xo O O O O O O O 6 o o qo O t 

jqoc ru! 1 0000 v^f><^ 0 (|oo 00006660000000000 ( t 5 tH> 0000006 d 0 ()OO 
^ O rrt I f h 0 0 0 op 0 0 6 O 00 O 0 000 Q O O 0 0 O 00 0 O Cl O 0006 D O O qo 6 1 

iCpcOni M n>OUdO6.i>p^666OOOOOOOOO0ODfH>OOr deoqoppqpQOunqo 
09000000000 0 UXIp OOOpO 0600000 000 O^CTO’OO'O 0 0 0 OX>d 0 00 0 O O do O i 
>c|oo()ooooooooopooo 6 ooooooooooo 6 >ooooo 66 oooooqoocKv>oooc 
O 6 0 0 o D 0 000 O O 00000000 00000 0 0 OOOt) 0600000 O O O 0 O 60 T) (>O 0 0 0 0 ( 


OUACHITA MOUNTAllNS 


ARKANSAS 


test 

Site 


10 0 


Kilometers 


Figuj^e 2 * IiOGStxon Of tjGst sit^e. 


GF'Vjtm‘1 »• ..C>, 

)Ci: ^ ii/r*, : 




SGgmeivta) which, on the west, are coincident in part with 
mapped normal faults (Fig. 3). 

The Ouachita Mountain fold belt, which is south of the 
Arkoina basin, is characterized by contorted imbricate 
thrust plates, intense folding with thrust faults, moderate 
folding without thrust faults, and gentle folding (Diggs, 

1961) . 

All rocks exposed within the test site are of early 
and middle Pennsylvanian age and are within the Atokan, 

Des Moinesian, and Morrowan Series. The Atoka Formation of 
the Atokan Series is essentially a sandstone, siltstone, 
and shale sequence and comprises 95 percent of the surface 
exposures within the test site. In the subsurface at least 
35 different gas-'bearing zones have been recognized, rangi* ^ . 
in age from early Pennsylvanian to late Ordovician; however, 
the principal gas-bearing zones are sandstone units within 
the Atoka Formation. 

Although much of the ear 3.y drilling and production 
(early 1900s) was associated with surface-expressed structural 
traps (anticlines) , more recent drilling and subsurface 
studies provide evidence thaf» l®est for the Atoka 
gas-bearing zones in the central part of the Arkansas part 
of the Arkoma basin, there is little discernible relationship 
between structure and gas production. Variation in porosity 
and permeability, i.e., stratigraphic entrapment , appears 
to be the primary factor influencing the presence of gas. 
Structure is generally of only secondary importance. 
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Figure 3. Landsat mosaic, winter scene showing prominent northeast- 
southwest-trending topographic lineaments. 







The Arkoma basin is essentially a dry gas province (no 
liquid hydrocarbons associated with gas production) , the 
gas being composed of approximately 95 percent methane by 
volume. At the end of 1978, there were 71 gas fields in 
the Arkansas part of the Arkoma basin. During 1978, 
89,373,839 MC)?G (1000 cubic feel- of gas) was produced from 
an estimated reserve of 1,542,739,839 MCFG (Arkansas Annual 
Oil and Gas Keport, 1978). 

111. PREVlOaS INVESTIGATIONS 
A. Linear Terminology 

Practures, fracture traces, microfractv;a?es,’ meso- 

fractures, macrofractures, megnjoints, lineations, linears, 

■ ■ ■■ 

and lineaments are some of the names given to naturally 
occurring, straight-line, aligning features observable by 
remote sensing. Many investigators have classified linear 
features according to their length (Blanchet, ’ 1957? Lattman, 
1958; Ilaman, 1961; Hoppin, 1974; Johnson, 1974) . These 
investigators found that the relatively long linear features 
are regionally significant whereas the relatively short 
features are more closely related to local conditions. 

Dlanchet (1957) describes fractures as ''generally 
abundant, natural lineations discernible on aerial 
photographs . " Fracture traces are defined by Lattman (1958) 
as natural linear features that have less than 1. 6 km 
continuous expression. Lattman terms fractures longer 
than 1 . 6 km "lineaments . " Haman (1961) uses the term 


"inicrofracture” to describe linear featnros less than 3,2 km 
long and the term ''macrofracture'' for features more than 
3,2 km long. Johnson (1974) describes regional lineaments 
as relatively long linear features which have been propagated 
upward from the basement rock. His boundary or local 
lineaments, in contrast, are relatively short and wore 

\i {/ 

generated by differential compaction along the flanks of a 
subsurface anomaly , Hoppin (1974) differentiates laetv^reen 
lincars and lineaments: -'Linears are single reetilinear 

elements commonly, but not necessarily, of sfcructviral origin. 
Lineaments are generally rectilinear lines or sones of 
structural discordance of regional (100 km or longer, ) 
extent , *' 

An excellent terminological reappraisal of lineaments 
and linears is provided by O'Leary et al, (1976) . These 
authors suggest that "lineament'' be used instead of 
"linear” for describing simple or composite linear features. 
Werner (1976) similarly sujnmarises the use and misuse of 
linear and lineament terminology. He suggests that a linear 
is similar to a fracture trace or lineament that is 
continuous , whereas a lineament can bo composed of an 
alignment of discrete features and short segments of lines. 

The term "linear" is used in this investigation 
because, by definition , it does not apply to composite 
features whose parts may bo aligned. A linear, as used in 
this report and inferred from remote sensor data, is a 
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natuiralXy occurring » continuous straight-line features, 
regardless of length, which may be manifested as a tonal, 
vegetative, stream pattern, or landform alignment. Two 
alig'nev but distinctly separate linears would be mapped as 
two entities (linears), •having the same azimuth (compass 
orientation) and perhaps different lengths. 

B. Linear Analysis and Geologic Significance 

Since the launch of Landsat-1 in 1972, considerable 
interest has been generated in linear feature tectonics, 
especially as related to regional fracture patterns. Studies 
conducted throughout the world have been reported in two 
International Conferences on New Basement Tectonics 
(Hodgson, 1976; and Podwysocki and Early, 1976). Geologic 
applications of fracture, linear, and lineament analyses 
derived from various remote sensing techniques have had 
varied success. 

Many workers involved in geologic remote sensing analyses 
would agree that global patterns Of regional linear features 
are related to pre*'CKis ting zones of weakness reflected into 
younger sedimentary strata by upward propagation, probably 
from the basement complex, whereas "local" linear features 
result from nonregional stress conditions indicative of 
subsurface discontinuities (Blanchet, 1957; Mollard, 1957; 
Hodgson, 1961a, b; Wobber, 1967; Wise, 1968, 1969; Rumsey, 

1971; Haman, 1972; Johnson, 1974; Kowalik and Gold, 1975; 


11 


Sanders and Hicks, 1976). Similarly, most structural remote 
sensing studies reveal that linear features are commonly 
associated with zones of increased joint or fracture 
concentration (Hobbs, 1911; Lattman, 1958; liattman and 
Matzki, 1961; l*attman and Parizek, 1964; Overbey et al., 
1974). The degree of success and geologic utility of remote 
sensor linear feature analysis, however is, a subject of 
considerable debate. Fracture trace, linear, and lineament 
analyses for groundwater exploration have been moderately 
successful in specific geologic environments. Significantly 
greater yields were obtained from water wells located in 
proximity to linear features mapped from remote sensor data 
(Lattman and Parizek, 1964; setzer, 1966; Sonderegger, 1970; 
Siddiqui and Parizek, 1971; Moore, 1976). However, not all 
hydrologic studies produced positive correlation between 
linear features and water yields (Vvhitesides , 1971; Ogden, 
1976) . 

Several articles have outlined the importance of 
recognizing linear features as a possible petroleum 
prospecting tool; however, the success or failure rates of 
actually using features identified by remote sensing 
techniques is an exploration program have not been well 
documented. Saunders and Hicks (1976) provide a summary 
of articles showing the importance of wrench faults in 
creating potential petroleum traps. These authors suggest 
the application of satellite lineaments in prospecting for 


concealed wi:onch faults oir other basement faults which may 
control massive reef growth, Attempts to predict petroleum- 
bearing rock fractures from surface linear feature analysis 
have had mixed success, Some authors (Overbey and Rough, 
1971; Aipay, 1973) suggest that linear feature analysis used 
in conjunction with joint strike measurements and borehole 
data can provide significant information on the orientation 
of natural fractures in a reservoir. The mapping of linear 
features has been shown to be helpful in exploration for 
some fracture reservoirs (Rumsey and Gelnett, 1976; Ryan, 
1976; Zirk and Lahoda, 1978), However, other auth°^‘® 
the contention that although linear feature analysis may 
serve at least partially to predict subsurface fracture 
orientation, the additional step of predicting a relationship 
between petroleum production and proximity to subsurf ace 
fractures is less certain (Overbey et al., 1974; Werner, • 
1978). Jones and Rauch (1978) and Beebe and Rauch (1979) 
give evidence that Iiandsat lineaments are associated with 
low yields from gas wells into Devonian shales, yet short 
photolinearaents have a positive correlation with high gas 
yields. Jackson et al. (1979) could find no basis for the 
validity of remote sensing in the exploration for shale gas 
in West Virginia. 

Fractures and faults in hydrocarbon-bearing strata have 
been recognized as factors affecting migration , accumulation, 
and production of hydrocarbons. Carbonate and less permeable 
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sandstone iresorvoirs are particulax'ly susceptible to 
fracturing. Efforts to determine the nature and extent of 
subsurface fracture systems and their effect on hydrocarbon 
reserves and reservoir flow during production have generally 
been unsuccessful. The complexity of establishing an 
acceptable relationship between surface and subsurface 
features is presumably duo in part to operator variability, 
different tectonic regimes in various study areas, differences 
in remote sensor data sets, non-standardization of analysis 
techniques, mai'ked scale diffeiencos in imagery data bases, 
and probably most important, lack of understanding of the 
possible origin of the mechanisms responsible for the 
observed features (Podwysocki and Go3.d, 1976). 

IV. SURFACE JOINT STUDIES 
A. Data Collection 

Joint data wore collected at 67 different stations 

2 . . . 
throughout the 11,000 km study area. In addition to joint 

orientation and frequency, bedrock lithology, thickness, 

and inclination were recorded at each station. At most of 

the stations the dip of the joints is essentially vertical. 

As illustrated in Figure 4 the regional orientation of 

surface joints within the study area is generally nearly 

north-south and east-west; however, deviation from the 

regional pattern can be seen at stations 19, 27, 28, 45, 

46, and 62. 
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The characterletio two sets of orthogonaXXy oriented 
joints found at each outcrop are iXXustrated in Figure 5. 

One set* composed of straight* through-»going joints* is 
designated the systematic set** its average orientation is 
NX0®W. The other set# which is more irreguXar and contains 
joints that coiwitonXy terminate against the systematic joints, 
is designated the non-systematic set; its average orientation 
is NSe^B. This terininoXogy is modified from that used in 
previousXy pubXished reports by Hodgson (X961a,b). PXumosc 
structures were observed on the surfaces of numerous 
systematic joints (Pig, 6); however* plumose structures wore 

i 

not observed on the non- systematic joint surfaces. Plumose 
structures are believed to be indicative of tensional stress 
conditions vGramberg* 1966). 

B. Fracture System Origins 

Many areas such as the Colorado Plateau* 'Vinta basin* 
and Piceance Creek basin are characterized by one or more 
sets of regional orthogonal fracture patterns (Stearns and 
Friedman, 1972). Several authors have discussed the possible 
origins of regional orthogonal joint or fracture systems ^^^^^^^^^ 
(Blanche t* 1957 ; Price* 1959; Hodgson* 1961; Grantberg* 1966). 
Stearns and Friedman (1972) summarize the various opinions. 

Although there is some variation in the trend of the 
axial traces of folds along the Ouachita Mountain fold belt 
(Pig. 3) * the general east-west trend of folding indicates 



Figure 5. Examples ot ortnogonai jointing in the 
Arkoma basin. 




a riorthetly compressive stress. If this hypothesis is 
correct, the systematic joint sets in the study area are 
oriented parallel with tho major Ouachita compression and 
the non-systematic joirt sets are oriented perpendicular 
to the compressive stress* I’his pattex-n conforms well with 
theoretical models for vertical fracturing presented by 
Price (1959) and Gramberg (1966) , According to tho studies 
of Price and Gramberg, the stress eonditions needed to form 
joints similar to those the Arkoma basin are: (1) vortical 

major prinGipal stress, (2) north-south intormediato px'incipal 
stress, and (3) east-west least principal stre^ss. 

The basic mechanisffi for joint formatiutt in the Arkoma 
basin is suggested to be tensional stress caused by bux'ial 
and later northward 'compressiGn atti'ibutable to the Ouachita 
orogeny. The theories of Price (1959) and Gramberg (1966) 
each require the strata containing the joints to have been 
buried relatively deeply during the Gompressive phase of tho 
Ouachita orogeny* According to Diggs (1961), major Ouachita 
compression was not active until late Pennsylvanian time. 

The thickness of sedimentary rock within the Atoka Formation 
(early Pennsylvanian) ranges upward to at least 6,000 m 
(Branan, 1968). In some areas whore erosion lias not yet 
taken place, the Atoka Formation is overlain by more than 
900 m of younger sodimontary rock. It therefore seems 
probable that, at the time of major Ouachita compression , 
the strata comprising much of the Arkoma basin were subjected 
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to a very great vertical loadr sufficient to qualify as a 
major principal stress. Obviously, under certain conditions 
of extreme Ouachita compression, the major principal stress 
may have been reoriented north-south. However, this situation 
would have been conducive to the formation of regional 
shear fractures trending north-northwest and north-northeast 
(Price, 1/59), • conse^jhently , as soon as shear fracturing 
occurred, it would have reduced the north-south stress and 
restored the major principal stress to vertical. 

The evidence indicates that the basic mechanism for 
joint formation in the ArJeoma basin was tensional stress 
caused by burial and later northward compression by the 
Ouachita orogeny. Increased pore pressure also may have 
helped to initiate fractures by decreasing the effective 
stress in all directions, which would have resulted in a 
north-south tension fracture perpendicular to the east-west. 

t 

least principal stress . In addition, subsequent uplift and 
erosion of the Arkoma basin may have generated vertical 
tensional fractures as a result of the "residual stresses" 
stored in the rock in the manner described by Price (1959) . 

The non-systematic set of joints may have formed because 
of changes in the principal stress orientations resulting 
from the formation of systematic fractures. However, another 
explanation for the non-systematic joints is documented by 
Hodgson (1961a), who ;states: "Any warping or bending of the 

strata subsequent to the formation of the systematic joint 
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should generate stresses and, consequently, fractures in 
the blocks between the joints. The common occurrence of a 
right-angle or nearly right-angle juncture between the traces 
of non-systeroatic and systematic joints indicates some 
control exerted by the systematic joint on the forces 
producing the non-systematic joint . " 

C. Joint; Parameter Analysis 

Initial measurements were taken from each joint station 
on nine parameters: 

(1) Thickness of the bed that contains the measured 
joints (THICK). 

(2) Dip of the beds (DIP). 

(3) Orientation of dip (DIPDIR) . 

(4) Litliology of the joint medium (LITH) . 

(5) Frequency of the systematic joint set per 30 m of 
traverse measured perpendicular to the systematic 
jointing (SFREQ) . 

^ ^ of the non-systematic joint set per 30 m 

of traverse measured perpendicular to the non- 
systematic jointing (NSFREQ) . 

(7) Orientation of the systematic joint set (SORIENT) . 

(8) Orientation of the non-systematic joint set 

(NSORIENT) . fr' ' 


(9) Relative stratigraphic position of each joint 

station outcrop (ROCKUNIT) . (Older or lower rock 
units were assigned lower numbers.) 

From four of the parameters, numbers 5, 6, 7, and 8, three 

other parameters are generated: 

(10) Overall averaged frequency of the systematic and 
non-sysmatic sets (AVEGREQ) . 
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(11) Ratio of the systematic frequency of the non- 
systematic frequency (FREQDIF) . 

(12) Angle between the systematic and non-systematic 1 

sets (ORIENDIP) . 

Each joint station was categorised into one of three | 

provinces within the study area (Fig. 7). From north to I 

south, they are; I 

S ' I 

Province l; Area of high relief along the Boston | 

Mountain escarpment--characterized by I 

growth faulting artd minor folding. | 

, II 

Province 2: Area of low relief between the Boston | 

Mountain escarpment and the Arkan^sas | 

River — characterized by growth faults | 

and folding, | 

■ 4 

Province 3j Area south of the Arkansas River — I 

characterized by intense folding and j 

some thrust faulting. ' 1 




A correlation coefficient and significance probability 
level (a) were generated for each combination of the 12 
parameters. A correlation was assumed to be present if a 
relationship had a significance probability level greater 
than 95 percent (a < .05). Previous investigations have 
shown a correlation between joint frequency and other 
parameters such as stress, lithology, or bed thickness 

i\ .... 

(Price, 19i,>9; Harris et al., 1960; Knoring, 1965; McQuillan, 
1973). Consequently, the relationships of major interest 
in this investigation are between the joint frequency parameters 
(SPREQ, NSFREQ, and AVEPREQ) and the rock parameters (THICK, 
LITH, or DIP) . DIP is interpreted as an indicator of stress 
levels in most instances. 

Correlations are found between joint frequency and two 
otlier parameters — dip and, more significantly, stratigraphic 
position. No correlation is found between joint frequency 

and lithology or bed thickness. The dip of the beds has 

11 

a positive correlation with the average joint frequency at 
an a level of 0.043. One can infer from this relationship 
that the non-systematic joint set has a better correlation 
(a = 0.024) with dip than does the systematic set (a = 0.087). 
This relationship supports Hodgson's (1961a) theory that 
non-systematic joints are created by the bending or warping 
of the strata after the formation of the systematic joints, 
which produces fractures that are usually perpendicular to 
the systematic set. 


According to the joint origin theory of Grambcrg (1966), 
great depth with high confining pressure is characterized 
by a high frequency of relatively short, regularly spaced 
joints. As confining pressure is reduced (i.e., at higher 
levels) , the joints becoino less frequency and longer, and 
conunonly are in multi~frachure zones. In this study the 
relative stratigraphic position of each joint station 
outcrop shows a negative correlation with joint frequency 
(« = 0.012). In other words, the joint frequency increases 
in the lower part of the stratigraphic succession. This 
finding corresponds to the result predicted by the Gramberg 
(1966) model . 

' *■ 

Another relationship of significance is the very good 
correlation (a = 0.0001) between joint frequencies of the 
systematic and non-systoniatic sets. This finding implies 
that the two sets are probably genetically related. A 
negative correlation (a » 0.006) between average joint 
frequency and province indicates that joint frequency decreases 
southward. This finding was not expected because the 
compressional stress component increases southward. However, 
joint stations in the southern part of the test site, on the 
average, are higher in the stratigraphiG succession than 
joint stations in the north, Therefore, the shift in joint 
frequency may be due to differencGs in the stratigraphic 
units measured, rather than the location of joint stations 
within the study 


In summary/ joints in. the Arkoraa basin consist of two 
orthogonally oriented sets of vertical joints. One sot 
(systematic)/ trending N10®W/ consists of straight/ planar/ 
through going joints. Its orienation perpendicular to the 
fold axes of the Ouachita Mountains on the south indicates 
that this joint set formed parallel with the major horizontal 
comprossivo stress. Plumose structures on many of these 
joints suggest a tensional origin. The other set of joints 
(non-systematic) has an average orientation of N86“E and 
consists of short and commonly curved joints which cut the 
rock between the systematic joints. The non-systomatic 
joints probably formed after systematic jointing in response 
to subseguent warping of the strata as indicated by a 
significant inGrease in non-systematic joint frogviency with 
increasing dip. 

V. linear REMOTE SENSOR ANALYSIS 

The end points of each linear from primary remote sensor 
data sets were digitized on a latitude (X) -“longitude (Y) 
coordinate system. Digitization was accomplished by means 
of a Tektronix model 4954 digitizing tablet which has a 
40-inch by 30-inch sensitive surface area and provides 
resolution of 0.01 inch. The tablet was attached to a 
Tektronix 4015 graphics terminal which/ in turn, was linked 
to an IBM 370-155 host computer. Data stored on magnetic 
tapes were used to produce linear and well location plots on 
a Zeta 3653 plotter. 
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A. Side-Looking Airborne Radar (SLAR) and Stereo Photography 

1. Existing Radar Data (SLAR) 

Parts of the study area had previously been imaged at 
various times by military APS-73 (x-band) and Westinghouse 
AN/APQ-97 (Ka-band) systems. Variations inflight direction 
(radar look-direction ) , altitude, and scale precluded the 
mosaicking of these flight lines. Similarly, because of 
marked differences in radar depression angles among the 
various flight lines and "spotty" coverage within the test 
site proper , linears obtained could be compared only in a 
visual , nonstatistical manner. 

2. Aerial Photography 

Stereo black and white photographs (generally 1:20,000 
scale) were obtained from various government, state , and 
private agencies. Approximately 1600 photographs were 
needed to provide stereo coverage . Unfortunately, the 
photography spanned some 20 years and the quality was not 
uniform. Mapped linears were transferred from the 

photographs to 7.5' orthophoto-quad bases and then digitized . 

2 

Because of the large area involved (11 , 000 km ) , 
different parts of the test site were assigned to different 
interpreters . In areas of interpreter overlap, considerable 
interpreter variability was apparent, especially in those 
areas in the low-relief Arkansas River valley. To provide 
an indication of the extent of interpreter variability, 
four interpreters were selected to analyze stereo photographs 
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in two diffoirent areas (thrG<3 7.5' quadrangle maps) oi; 
contrasting terrain configuration. The terrain of the 
Mountainburg and Mountainburg S.W. quadrangles is character- 
ized by high relief (>200 m) because of dissection of the 
nearly flat-lying strata. The Alma quadrangle , in contrast, 
is characterized by much less relief (<50 m) and includes 
part of the Arkansas River floodplain. 

Rosette diagrams for each quadrangle and operator 
(Pig. 8) indicate the number of linear s mapped (n), the 
average length of each linear (x) , and the percentage of 
linears mapped by more than one operator. The rosette 
diagrams drawn by each operatbJ-" show 'considerable variation. 
However, in the Mountainburg and Mountainburg S.W. quadrangles , 
the variation was less than in the Alma quadrangle. For 
example , in the Mountainburg quad all four of the operators 
identified a major linear set trending N42°W anv,, three of 
the operators approximated two other linear sets trending 
N23®E and N54®E. There is somewhat less agreement among 
rosettes constructed from the Mountainburg S.W. quad where 
three operators agreed on two linear sets trending N5°E and 
N25®E. Very few linears wore mapped in the Alma quadrangle 
and the linear trends recognized by different operators show 
no agreement. The number of linears mapped by all operators 
is by far the lowest in the Alma quadrangle. This finding 
may indicate that operators are more likely to distinguish 
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Figure 8. Results of operator variability study — stereo 
photography. 
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topographic linears than tonal or vegefcation-rela ted 
linear 5 caused by soil inoisturG variations. 

B. Radar (gAR) and RnhanGed Landsat 
1. Coinputcr l?rooessing and Enhancement 

In support of this rfeseracb effort, the Jet Propulsion 
Laboratory provided SAR iinagery of the test site obtained 
during May of 1977 and April of 1978 , Because of funding 
limitations, only the 1978 data set was digitally processed. 
This L*^band (2 3 -cm wavelength) imagery was obtained in a 
series of north-south flight linos with approKimately 8 -km 
_ swath widths and providing dual’-polari!&atiQn fioveraga. Two 
different Landsat scenos (winter and summer) were selected 
to provide maximum changes in solar illumination (elevation) 
and terrain cover. Both radar and Landsat imagery was 
subjected to various computer enhanGement and analysis 
procedures performed at the Imago Procassing Laboratory of 
the Jet Propulsion Laboratory. 

a. Radar Image Processing 

Strips of radar imagery which had been converted from raw 
data to image format by the JPL optical processing systeni 
were digitized on the IPL microdensitometer. Samples of 
film density were taken With a 50-micron aperture at intervals 
of 50 microns, corresponding to dimensions of 15. and 25 m 
in the range and azimuth directions r respcctivqlv. The 
resulting density measurements were scaled to teger 


values, although in any strip considerably fewer than this 
number of levels were usually occupied. A typical strip 
of sampled imagery yielded a digital array of apprOK’; 'tely 
500 samples in the range direction by 3500 lines in *» • 
azimuth direction. 

Removal of gain variation 

The radar system which generated the imagery produces a 
perceptible gain variation with range. This effect is the 
result of three phenomena : the normal spreading loss of power 

density vi^ith increased path length , the decrease in differ^ 
.ential bacjcscattering cross section (p**) as incidence angle 
increases, and antenna gain variations over the imaged swath. 
To combat this problem the system incorporates the capability 
to increase the gain as a function of range. The system is 
only partially effective in compensating for these effects, 
leaving a noticeable average intensity variation with range 
in the data. Such a residual effect is particularly apparent 
in a mosaic where the near range of one strip is adjacent to 
the far range of its neighboring strip. There is an obvious 
average intensity discontinuity at the strip boundary. 

To reduce the intensity descrepancy between adjacent 
strips, a scene-dependent gain correction was applied to 
each strip. The procedure consisted of applying a gain and 
offset to each line of constant range in the imagery to force 
the mean and variation (i.e. , ADS Ix^ x] where x^ is an 
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intake sample value, and x is the average o£ over a line 
of constant range) at each range to predetermined constant 
values. This process is admittedly inferior to a true cal- 
ibration operation, but the information available was in- 
adequate for a meaningful calibration. A major drav/back 
of the technique is the fact that the applied gain and 
offset might reflect the feature content of the imagery 
rather than correct for true gain variations. To reduce 
undesirable feature-cued gain modification, the process 
was controlled by using statistics obtained only from 
uniform-appearing areas in the mountainous regions. This 
method was generally successful in reducing apparent 
intensity jumps between strips in the regions where the 
statistical estimates were obtained. However, intensity 
discontinuities did appear in areas of low return where 
the gain attenuation was limited by the low signal strength. 
Because these areas were not of primary interest in the 
investigations, the artifacts were deemed acceptable. 

Conversion from range to surface distance 

Because the side- looking radar system records the 
returned signal from each pulse as a function of time, the 
imagery obtained from the raw data displays distance along 
the flight path in one coordinate (azimuth direction) but 
line of sight distanGe to the surface in the perpendicular 
coordinate (range direction) . A difference in range between 
two points nearly beneath the sensor corresponds to a much 


cjjroate)? sopairation on the suirfaco Ubon does an equal range 
clifference between points located to the side o£ the flight 
path. As a result, features appear greatly compressed in 
the near range areas of the imagery. 

This basic distortion is complicated by the fact that 
tho range to a surface point is dependent upon its elevation, 
V-Jitbout exact knowledge of the topography and the position 
of the sensor, an accurate conversion to distance on the 
surface cannot be made . The near range compression effect, 
however, is sufficiently pronounced that application of an 
approximate "first-order'' correction is often warranted. 

Such v^as the case in this project where the radar strips 
were to bo incorporated in a mosaic and registerod to Lnndsat 
imagery. The presence of gross local distortions would 
confuse the analyses and significantly degrade registration 
accuracy. Therefore, each strip of imagery was subjected to 
a geometric transformation to accomplish an approximate 
conversion from slant range to distance on the surf ace. The 
altitude was assumed constant over the strip. A simple trig- 
onometric conversion from measurement along the hypotenuse of 
a right triangle to distance on the horizontal side v;as made 
with the assumption that the topography was absolutely flat. 

Mosaic generation 

Radar strips which had been processed to remove gain 
variations and to convert from measurement of range to 


distoncG on the oui'faae were ready for incorporation into n 
mosaic, A mosaic was constructed by combining neighboring 
image pairs into two-component mosaics. These two-component 
mosaics were then merged in pairs to form four-component 

mosaics and the process was repeated to generate the 

* 

completc'd mosaic. At each step a pair of images were 
combined by deforming tlie right component to spatially match 
the left component at the boundary botwoen the tv?o images. 
Identical features were located manually in bands across the 
overlapping areas of both images . Typically, five bands 
of such “tiepoints" were obtained at equally spaced intervalB 
down the strip of overlap. To avoid compounding the 
complojdty of subseguent mosaicking steps, the right-side 
component was subjected to a strictly controlled geometric 
transformation to achieve alignment with its left-side 
counterpart. The trahSformation v/as performed in such a 
manner that each line of constant azimuth (a horizontal rov; 
in a component strip) was translated only in the horizontal 
direction and was translated as well as rescaled in the 
vertical direction. Between the bands of tiepoints the 
deformations were computed by linear interpolation between 
those indicated by direct observation at the tiepoints. The 
two-clement mosaics were formed by truncating the left and 
registered right-side components at the boundary and abutting 
the resulting strips. 
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Birightness transformations 

To provide image displays in which features ot^ various 
textures and dimensions could be easily discerned, the 
mosaics v/ere subjected to two brightness transformations. 
Providing the lower contrast was a transformation to produce 
an approximately Gaussian distribution of intensity. A 
transformation table was constructed to change the original 
brightness levels to the set of new values which produced 
the closest approximation to a Gaussian distribution with 
mean at midrange and six standard deviations spanning the 
256 available shades of gray. 

A more dramatic contrast enhancement was achieved by 
the second procedure. As in the first method, a transfor- 
mation table was generated to produce a specified distribution 
in the transformed image. In this case that distribution was 
uniform. The algorithm operated by the equilvalent method 
of forcing the cumulative distribution function to be linear. 

■ s.,:/';. 

spatial filtering 

The nature of the investigations conducted in this 
project required the clear visibility of fine detail in the 
imagery. To increase the apparent sharpness of the pictures , 
a spatial filtering operation was performed. This procedure 
increased the contrast of small scale brightness changes 
and reduced brightness differences between large scale 
features. If the low frequency (large scale) components had 
been completely removed, the filter would have been a 
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''high-p£vss'' liitej;, In this project contrast was 80 
percent attenuated in features larger than approximately 
2500 m and smaller dimension variations were loft unchanged, 
This contrast corresponded to bho average of a "high^pass*' 
filtered version with 20 percent of the original image. 

RGgistration to Ijandsat 

Any combined analysis of radar and landsat imagery is 
facilitated if a specific feature appears at the idGntical 
coordinatos in each typo of image, such spatial alignment 
(registration ) also makes possible the goneration of 
numerical eomposites of the images. To achieve rGgistration 
in this project the radar scene was geometrically manipulated 
to coincide with the Landsat scene. The transformation was 
controlled by approximately 80 tiepoints manually obtained 
by locating identical featurea in both images. A polynomial 
goometric transformation was used to describe the distortion. 
It was "fit'' to the tiepoints by adjusting its coefficients 
to minimise the sum of the squared error in its predictions 
of thetiepoint locations. To speed computation the polynomial 
transformation was applied only eit a sparse grid of locations 
in the image. Movements of points falling between the grid 
intorsoctions were obtainod through bilinear interpolation 
between those of the four surrounding grid locations. 


b. Land sat Image Processing 

Landsat imagery was obtained from the U.S. Geological 
Survey in standard digital form. Each scene consisted of 


images in four spectral bands sampled with a spacing of 


approximately 60 m across the scene and 80 m down the scene. 
Saroplcd at these intervals, the roughly 190 km square area 
imaged in a scene comprised a digital array of about 2300 
lines each containing approximately 3200 samples, 

Geometric correction 


The standard format of the USGS-supplied digital data 


is incompatible with much of the software in the ‘Image 


Processing Ijaboratory . Therefore , it was necessary to 
reformat the data as the initial processing step. During 
this "logging" process it was convenient and computationally 
efficient to undertake a set of geometrio corrections which 
could be performed independently on each image line. Included 


were the following geometric manipulations v 

(a) correction for the earth's rotation during the 
^ ^ ^ imaging the first and last 

linos , 


(b) 

(c) 

(d) 

(e) 


resampling to achieve equal horizontal and 
vertical sampling spacing, 


r enidval of synthetic samples introduced by NASA 
to ensure equal line lengths, 


correction for non-uniform sensor scan mirror 
velocity. 


correction for panorama effect, and 


(f) adjustment of sample coordinates to compensate for 
variation in alignment of the sensors of each 
spectral band. 
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Brightness ti:ansformatic^ns 

/ ' 

As in radar imagery# the distribution oJ: brightness 

// 

in the original Landsat images d|id not necessarily provide 

O 

sufficiGnt contrast and average level to produce useful 
photographic representations. To generate more informative 
visual displays# the individual spectral band images were 
subjected to the same intensity transformations that were 
applied to the radax- imagery# producing Oaussian and uniform 
distributions of brightness. 

Spatial filtering 

Edges and fine structure in the Landsat scenes were 
enhanced by tlie same filtering procedures that were applied 
to the radar imagery. v/ ' 

; / 

Color composites Ar 

Two types of color composites# referred to as standard 
and enhanced versions# were generated from Landsat imagery. 
The standaard version was obtained by modulahiag each of the 
throe primary colors (red# green and blue) with a single 
spectral band image. To ensure that spectral variations in 
all the bands v^ere equally displayed, the images were first 
subjected to the Gaussian brightness transformation . This 
process produced equal first- and second-order statistics i 
in the color component images and elimina ted the po s s ib ili ty 
that a single hue would dominate the scone . 


The enhanced color coniposiUe was produced to display 
even small spectral variations as marked color differences. 

The color enhancoraent procedure was designed to spread the 
distribution of brightness, causing it to occupy a broad 
region in the color sxjacci This effect was acconiplished by 
rotation to the space of oigenveGtOrs of the original R, G, 

B (rod, green, blue) distribution, Statistically indo.pendent 

images were forwod by assigning brightness based 

coordinates along axes aligned with each of the eigenvector s. 
Therefore, when a Gaussian brightness transformation was 
subsequently applied to each of the new component images, 
a roughly Spherical distribution was produced, generating 
a broader range of colors than the cigar-shaped original 
distribution. To maintain general color similarity to the 
standard color composite, the inverse rotation bac)c to the 
original color space was performed to generate the tri-stimulus 
components of the color image. 

Registration of multiple scenes 

To facilitate temporal comparisons of Land sat imagery, 
scenes obtained at different times were registered. The 
procedure was identical to that used to register radar to 
Landsat imagery. The accuracy of the LandSQ^ 
however, was considerably greater because the geometric 
differences were less severe and more continuous than the 
rather disjoint distortions generated by mosaicking 
independently correct radar strips. 


c. Radar /Land sat Composite 

... m 

Registration of a radar mosaic to a Landsat scene made 
possible point-by-point manual comparisons of the two data 
types. It also allowed numerical combinations to be 
undertaken. For this project composites were produced by 
using Landsat data to provide color information and radar 
imagery to supply brightness modulation. The enhanced 
color tri-stimulus components were transformed to a spherical 
coordinate system known as hue saturation and intensity 
(HSI) space. The intensity coordinate is a distance along 
a vector forming equal angles with the orthogonal R# G, B 
(red, green, blue) axes. Hue is the angular rotation in a 
plane perpendicular to the intensity vector. Saturation Cs 
the angular displacement at the R, G, B origin away from the 
intensity vector* In this system brightness information 
is contained in the intensity coordinate and color in the 
hue and saturation dimensions. To form the composite, the 
Landsat intensity image was simply replaced by the registered 
radar mosaic. The inverse transformation from HSI to RGB 
space was then performed to obtain the tri-stimulus 
components of the color composite. 

2. Evaluation of Enhancement Techniques 
a. Landsat-Band Comparison 

The first evaluation is for the four bands of each 
enhancement and scene combination. For this evaluation, 
total number of linears detected by each of three different 


operators (A, B, C) is used. The Landsat sceno-enhanceinont 
combinations are: 

Landsat winter - Gaussian stretch 

Landsat winter ~ uniform distribution stretch 

Landsat summer ~ Gaussian stretch 

Landsat summer - uniform distribution stretch 

Tables 1, 2, 3, and 4 show the number of linears drawn 
by each operator (A^ B, C) for each band and the band and 
histogram numbers (Appendix) . The highest nuinber for 
each operator is indicated by * and the lowest number by x. 

Table 1. Landsat winter — Gaussian stretch. 



(1) 

Band 4 

(2) 

Band 5 

(3) 

Band 6 

(4) 

Band 7 

A 

127 

178 

101 

229* 

B 

108 

212* 

84^ 

103 

C 

33 

28^ 

62 

63* 

Table 

2. Landsat 

winter — unifojrra distribution stretch. 


(5) 

(6) 

(7) 

(8) 


Band 4 

Band 5 

Band 6 

Band 7 

A 

133 

155* 

119^ 

127 

B 

114 

209 

113^ 

269* 

C 

47^ 

54 : ^ ^ 

53 

75* 


Table 3. 

Landsat summer- -Gaussian 

stretch . 


(9) 

Band 4 

(10) 
Band 5 

(11) 
Band 6 

(12) 
Band 7 

A 

74 

59X 

86* 

86* 

B 

130 

’ 147 

95^’ 

102 

C 

16 

15^ 

62* 

40 

Table 4. 

Landsat summer — un i f or m 

distribution stretch. 


(13) 

(14) 

(15) 

(16) 

Band 4 

Band 5 

Band 6 

Band 7 

A 

58^ 

124* 

100 

93 

B 

143* 

126 

134 

124^ 

c 

21 

20^ 

30 

51* 

For 

the Landsat 

i winter scene, 

with both (Saussian and 


uniforru distribution stretchy two of the three operators 
found the greatest number of linears on band 7. For the 
Lands at summer~*-Gaussian stretch combination ^ two of the 
three operators found the greatest nuitiber of linears on band 
6. Hone of the operators agreed on which band had the, 
greatest number of linears for the Landsat summer — uniform 
distribution stretch combination. 

For the Landsat winter scene, regardless of enhancement 
technique, two of three operators found the greatest number 
of linears on band 7 . There was less agreeroent among the 


opera tojrs £ot the Lanclsat summer scene. In one case, band 
5 was the best and in the other case there was no agreement. 

In two of four cases the least number of linears was 
found on band 6. However, in three of four cases, band 6 
showed the least Variation in the number of linears drawn 
by the three operators. 

Also, from Tables 1, 2, 3, and 4 it is apparent that 
the Landsat winter scene detects more linears than the 
Land sat summer scene. The uniform distribution stretch 
consistently detects more linears than the Gaussian stretch. 
The best combination to find the gi'eatest number of linears 

r 

is Landsat winter, band 7, uniform distribution stretch. 

b. Enhanced Color — Standard Color Comparison 

For this evaluation, the number 'of linears detected with 
enhanced color Landsat imagery is compared with the number 
of linears detected on standard color Landsat* imagery. The 
Landsat scene-enhancement combinations are: 

Landsat winter — enhanced color 
Landsat winter — standard color 
Landsat summer — enhanced color 
Landsat s ummor — standard color 

Table 5 and 6 show the number of linears drawn by each 
operator (A, B, CO for each scene and the histogram number 
(Appendix) . The highest number for each operator is 
indicated by *. 
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Table 5, Landsat winter. 



(17) 

Bnhanced Color 

(18) 

Standard Color 

A 

178 

203* 

D 

191 

192* 

C 

51 

62* 

Table 6 . 

Landsat summer. 

p-'»^iiiiMriiMii lim vmn'ian* ofciwwiriiiiyiiiwiiii .i,..,. 



(19) 

Enhanced Color 

(20) 

Standard Color 

A 

143* 

129 

B 

148 

167* 

C 

37 

41* 


The standard color Landsat imagery detects more linears 
for both the summer and winter scenes. Also, a greater 
number of lineras is detected with the winter scene. 

C. Synthetic Aperture Radar — L-band 

The images from which histograms 21 and 22 were produced 
are uncontrolled mosaics whereas the images from which 
histograms 23 and 27 were produced are controlled mosaics 
(response range dependence removed) . Histograms 21 and 22 
wore obtained from the 1977 data set, Table 7 compares the 
number of linears from each of the three operators (A, B, C) 
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with tho histograms indicated (Appendix). The greatest 
number of Xinears found by each operator is shown by * and 
the lowest number by x. 

Table 7. Linear'-operator variation. 



(21) 

(nil) 

(22) 

(HV) 

(23). 

(IIU)'^ 

(27), 

(V)I)-^ 

A 

65 

50^ 

85 

111* 

B 

54 

49^ 

137* 

87 

C 

23 

4^ 

41 

100* 


^Response range dependence removed. 


Table 7 shov^;s that two of the three operators found the 
greatest number of linear s on the (VH) (response range 
dependence removed) image and all three operators found the 
least number of linears on the (HV) uncontrolled image. More 
linears are detected by mosaicked (response range dependence 
removed) April 1978 imagery than uncontrolled imagery (May 
1977). 

Enhancement comparison 

Three SAR (L-band) enhancements — uniform distribution 
stretch, high-pass filtering, and high-pass filtering with 
20 percent added back are evaluated. As before, the evaluation 
is based on the number of linoras detected. 

Table 8 compares the number of linears from each of the 
three operators (A, B, C) with the histograms indicated 
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(Appen<3ix) . The greatesit nunibar of Xinears found by each 
operator ie shown by 

Table B. Enhanced radar. 



(24) 

Uniform 

Distribution stretch 

(25) 

High Pass 
Filtering 

(26) 

High Pass 
201 Added Rack 

A 

60 

109* 

80 

B 

143 

72 

196* 

C 

82* 

80 

79 


Although none of the operators agreed on which type of 
enhancement detocts the most linoars, the high-pass 
filtering shows a very small variation in number of linear s 
found by the three operators. 

d . Radar /Landsat Merge • 

The parallGl-polarised (HH) {response range dependence 
removed) imago was registerod to the winter and suntmor 
Landsat scenes. Table 9 compares the number of linears from 
each of the three operators (A, B, C) with the histograms 
indicated (Appendix) . The greatest number of linears found 
by each operator is shown by * . 

All three operators found more linears on the radar/ 
Landsat winter scene than on the radar/Landsat summer scene . 
Unfortunately, because of the size limitation of the color 
playback device , the area generated for the radar/Landsat 
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TabXo 9, Radai* plus Landsat. 


(28) 

Win tor Scono 


(29) 

Summor Scono 


A 69* 40 

B 60* 46 

C 74* 33 


analysis was smallej; than that used for other linear analysas 
and histogram generation. Table 10 indicates modified area 
histograms for comparable areas derived from the. radar/ 
Landsat input combinations that provided histograms 28 and 
29. The greatest number of linears found by each operator 


is shown by*. 

Table 10. Modified area histograms. 

t tai" " I '! rt?7'rT i n'; i i t.; ^ I m j ;'fe';Trr:7'vwj w k r'li 

(32) 

(30) (31) KJI-Response 

Landsat Summer Landsat Winter Range Depondenco 
Enhanced Color Enhanced Color Removed 

A 95 162* 49 

B 79 138* 130 

C 27 42* 30 


When histograms 28 and 29 (v/inter and summer radar/ 
Landsat combinations) are compared with the original data 
(histograms 30-32 ) , it is obvious that the Landsat winter 
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enhanced scone provided for the detection of most iinears. 

The radar /Land sat merge appears to improve tho confidence 
of operators in interpreting relatively longer lincars. 

The average length of linears on the radar/Landsau merge 
shows a moderate incroase for two of the three operators, 
but the total number of linears is smaller, Gperator c 
was able to detect more linears using the radar/liandsat 
merge than with either the radar or Landsat alone. However, 
the other two operators show marked disagreement. Again » 
operator variability was considerable. 

Although the VIJ mosaic provided for the detection of 
relatively more linears, none of the SAiR enhancement 
techniques is significantly better than tho others, 

e. Radar and Radar /Landsat Evaluation 

From the preceding diseussion it is apparent that L-band 
radar imagery provided for detection of fewer lincars than 
the optimum data set, that is, tho Landsat winter scene. 

This conclusion is not surprising , however , if one examines 
the contrast in terrain displays provided by the two remote 
sensor data sets. The dominant linear features inferred 
by all interpreters on both tho Landsat and radar imagery 
were topographic linears , i . e. , usually alignment landf orms 
or straight valley segments. Topographic lincars arc 
enhanced on the Landsat winter scone because of the relatively 
low solar illuminatio;< (21® solar elevation) which causes 
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highlicjhting and shadowing on many o£ the terrain elements 
oriented toward and away from the sou theai.\ tern solar azimuth 

|i 

(142®, measured clockwise from north). This low sun angle 

enhancement is especially evident if the 'Winte*^ scene 

// 

(Pig. 9) is compared with a high sun angie (56® solar 

* I/* 

elevation) summer scene (Pig. 10). Thd| contrasting terrain 
information in the two Landsat scenes is striking, but not 
a new observation, in recent years many bandsat invostigatois 
have found that the vi?inter months' lower sun angles enhance 
topographic detail. However, the similarity between low 
solar illumination and radar illumination is most important 
in analyzing the ShR imagery linear detection. 

The L-band radar imagery was obtained along north-south 
flight lines, with two separate but overlapping swaths 
providing flight line coverage (Pig. 11). This flight line 
configuration allowed for a continuous change in incidence 
angles from 43° (near range) to 71° (far range) , or 47° 
to 19® depression angle, respectively. In contrast, the 
Landsat winter scene provided a constant solar elevation 
(equates to radar depression angle) of 21® across the entire 
swath. Only in the extreme far range part of each radar 
swath does the radar depression angle approximate or equal 
the solar elevation of the Landsat winter scene. The 
enhancement of topographic linears on the radar imagery is 
a function of two parameters, one related to the imaging 
system and the other related to terrain. The system parameters 
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are radar depression angle and look-direction. Changing 
the depression angle can have a significant effect on the 
enhancement of topographic linears. The May 1977 li-band 
imagery was obtained using a swath configuration that 
provided for much larger (steeper) depression angles across 
the swath (Pig. 12), Although the imagery was not processed 
in the same way as that of 1978, and was not analyzed in a 
controlled mosaic format, histograms 21 and 22 (Appondixf 
1977 data) show significantly fewer linears detected with 
1977 than with 1978 data (histograms 23, 24, 25, 26, 27). 

Other researchers have recognized this relationship 
(MacDonald et al, , 1969; Gelnett, 1977). Gelnett (1977) 
provides a summary of articles emphasizing the importance 
of considering radar depression angle and look-direction 
when attempting linear analysis from radar surveys. It should 
be noted, however, that because of east and west look-directions 
of the radar flight lines, some linears oriented parallel or 
nearly parallel with the flight lines were more easily 
detectable on the radar than on the Landsat imagery. Thus, 
in a way similar to radar, Landsat is also look-direction 
dependent (single solar azimuth illumination angle ) , Gelnett 
(1977) reports that single-look radars provide only about 70 
percent of the available geologic data. Therefore , because 
of radar ' s operational capability of selectable look-angles 
and depression angles, radar surveys should complement any 
topographic linear survey using Landsat. 






From the complementary aspects of Landsat and radar^ 
intuitively one might expect the radar/Landsat merge to 
provide better linear mapping capability than a single-sensor 
mapping mode. Comparison of histograms 28 and 29 (winter 
and summer radar/Landsat merge, respectively) with histograms 
30, 31, and 32 (original imagery used for Landsat-radar 
merge) provides no evidence of an overall improved linear 
detection capability for the merged data sets. Although each 
operator was able to detect linears of greater average length 
than those detected in the original data sets, the merge of 
the two data sets tends to "masik" the relatively smaller and 
more subtle linears. The enhanced color of the Landsat, 
for example, suppresses subtle radar linears that are 
relatively distinctive by demarcation between black and 
white tonal contrasts on the original imagery. A black and 
white radar/Landsat merge of Landsat uniform distribution ■ 
Stretch, band-7, winter scene with the VH radar appears to 
be an ideal data set. Such a black and white combination 
was not generated for this study. 

f. Operator Variability 

The problem of operator variability in the mapping of 
linears and lineaments has been addressed recently by 
Podwysocki et al. (1975), Siegal (1977), and Huntington and 
Raiche (1978). All found significant discrepancies in both 
coincidence and directions of linears drawn by different 
operators on the same Landsat image. Podwysocki et al. 


report that of all linear s mapped among four operators ^ 
only 0.4 percent were recognized by all four operators, 

4.7 percent by three, 17.8 percent by two, and 77 percent 
by one operator. Siegal <1977), using five operators and 
a single Landsat scene, found that only 4 percent of the 
lineaments were consistently delineated by all five operators 
Huntington and Raiche (1970) suggest that differences between 
two interpretations, by two observers looking at the same 
scene, can be of the same order as would be expected for one 
person interpreting two scenes with entirely different 
geology. 

Differences between linear interpretations of the same 
scene can be extraordinary. Even when the operator variables 
are minimized, contrasts between remote sensor data sets 
may produce significant differences in interpretations. 
However, maps derived from linear density analysis are found 
to be somewhat less influenced by operator variability. 
Generally, conservative linear "mappers" tend to infer fewer 
linears per unit area, but relative concentrations are 
sometimes agreed upon among operators. 

VI.. PROXIMITY ANALYSIS— LINEARS AND GAS YIELDS 

Production or gas yield data (initial calculated open 
flow) combined with well locations and linears from various 
remote sensor data sets were called up simultaneously from 
the Zeta 3653 plotter. Figures 13 and 14 are representative 


PRODUCTION PLOT 


LEGEND 

+ = NON-PRODUCING WELL 
A = < 2 MILLION CUBIC FEET 
^ “ 2 < 5 MILLION CUBIC FEET 

® = 5 < 15 million cubic FEET 
0=15 MILLION CUBIC FEET OR MORE 

NUMBER OF PRODUCING WELLS = 1059 


Figure 13. Linears from Landsat winter scene, band 7, 
IPCOF (initial production calculated open 
. flow) indicated. 
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PRODUCTION PLOT 
LEGEND 

+ • NON-PRODUCING WELL 
A = < 2 MILLION CUBIC FEET 

^ = 2 < 5 MILLION CUBIC FEET 
O = 5 < 15 MILLION CUBIC FEET 
° = 15 MILLION CUBIC FEET OR MORE 

NUMBER OF PRODUCING WELLS = 1059 


Figure 14. Lincars from SAR L-band radar and well 

locations with IPCOF (initial production 
calculated open flow) indicated. 
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sample plots fcom the Landsat band 7 winter scene and L-band 

respectively. The superimposed plots provide an ovotall 
coroparison of poor, average, and good gas well yields with 
respect to linear position, clustering, and density. 

Regardless of the remote. sensor data set used, it was apparent 
that the test site could be subdivided into two contrasting 
regions: the southern two-thirds of the test site which is 

characterized by an overall lack of linears in an area where 
producing wells are most numerous, and the northern third 
of the test site which has a relative abundance of linears 
where gas wells are nearly absent. An attempt is made to 
define the relationship between gas production and linear 
proximity in both of these regions. 

A. Southern Region 

Within the southern region or Arkoma basin proper, where 
wells are numerous and linears sparse, all attempts (linear 
regression analysis) to establish a relationship between 
well productivity and linear proximity proved futile. 

However, two additional methods of analysis were tried With 
well stimulation (hydraulic fracturing) data. Hydraulic 
fracturing, first introduced in 1949, is a process of 
injecting fluid into a wellbore under sufficient pressure 
to induce a vertical or horizontal fracture in the well 
formation. Continued injection of fluid after formation 
"breakdown" has occurred extends the fracture laterally from 
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the wellbore. Eventually an eguillbrium is reached in which 
the filtration fluid losses through the walls of the fracture 
are equal to the fluid injection rate, and lateral extension 
of the fracture ceases. It is generally assumed that if 
the fracture area so generated is to remain open to flow a 
proppant must be simultaneously injected into the fracture 
area. The most common propping agent is 20-40 mesh sand 
(0.42-0.84 mm), though other propping agents are used such 
as glass beads, plastic particles, and ground walnut hulls. 
Various completion data (amount of proppant, proppant size, 
and fluid injection volume) from well completion * cards were 
used to evaluate hydraulic fracturing' effectiveness. This 
effectiveness was then compared with linear proximity or 
linear density in a manner somewhat similar to that reported in 

shale-gas studies where surface joints or fractures were 

■//' 

used as predictors of the orientajbion of induced fracturing 
in the subsurface (Overbey et al., 1974; Jones and Rauch, 

1976) . Initially, a computerized dat^ system was generated 
having the following information taJcen from completion cards 
for each of the producing wells in the study area. 

Producing Formation (up to 3 
dif f erent zones ) 

Open Flow Potential 
Perforation Interval 
Acidizing Volume 
Sand Proppant Weight 
Hydraulic Fracturing Volume 
After Fracture Test Production • 

(not always available) 


Well Name 
Field Name 
Legal Description 
Longitude/Latitude 
County 
Spud Date 
Completion Date 
Total Depth 


different methods *of anlaysis were applied with well 
stimulation data, one method for linear density and the 
second for individual linear proximity. The first method 
was an attempt to ''generate*' additional linears by use of 
the Ronchi grating in the area of maximum well density. The 
Ronchi grating or square wave transmission grating consists 
of a set of opaque parallel lines (normally 200/inch) 
separated by equal-v/idth spaces on a transparent backing 
such as glass or plastic. Light is diffracted by the Ronchi 
grating in a direction perpendicular to the direction of 
the lines on the grid. When an image is viewed through the 
Ronchi grating the diffracted light will be spread out in 
one direction, but undisturbed in the perpendicular direction. 
This diffraction tends to blur out linears on the image that 
are parallel with the lines of the grid, thus accentuating 
the linears that are perpendicular to the lines of the grid 
(Pohn, 1970 5 Bench et al., 1977; Jackson, 1979). 

It Density Analysis 

■.r- 

A contour map of linear density based on the Landsat-2 
(winter scene) imagery was prepared. No significant trend 
or correlation could be established from this analysis , 

The percentage distribution of fractured and non-fractured 
wells was almost identical for each contour interval (see 
Table 11) . Also, no trend could be established with respect 
to the volume of fracturing fluid injected and the linear 
density. 


Table 11. Induced hydraulic fracturing linear density 
comparison. 


liinear Density 
(m/sq km) 

Fractured 

Wells 

Non-Fractured 

Wells 

No. of Wells 

% 

No. of Wells 

% 

0-176 

318 

(i 

73 

452 

72 

176-353 

91 

21 

119 

19 

353-529 

17 

4 

30 

5 

529-706 

3 

1 

8 

1 

706-882 

2 

- 

9 

2 

882-1059 

3 

1 

7 

1 


434 

100% 

625 > 

100% 


2, Proxiinity Analysis 

The second method provided the opportunity to call up 
wells in a given localized "window" enclosing a particular 
linear field or section of a field. The latitude and 
longitude of the two terminal points of a linear were read 
into the computer to create a rectangular window based on 
specified distances from these end points (Fig. 15). All 
wells within this window were plotted with an accompanying 
printout of the data parameters shown in Figure 15. Wells al- 
so could be called up by any one or combination of the five ma 
j or geological zones known to yield gas in the areas selected. 

Four separate areas were examined in detail by this 
method in terms of contrasting linear densities/ proximity 
of gas wells to major linears, and contrasting linear 
orientations. 
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DATA INPUT: Longltudc/latltude of linear end points, 

Dj^, D 2 , and scale choice, formation key, 

DATA OUTPUT: Plot shown above with well record printout 

as shown below. (D - denotes dry hole, other 
letters denote coding for producing wells,) 



WELL NAME = WILSON PROD GO. //I SYDNEY JOHNSON FLD NAME = DOVER j 

LEGAL DESCRPT = SEC 2-9N“20W,132O'FNL,26OO’FEL 

SPUD DATE= 02-01-74 LONG* 93DEG 5Min 24SEC LAT= 35DEG 27MIN 45SEC 

COMP DATE® 03-01-74 TOTAL DEPTH® 3971’ COUNTY® POPE | 

GAS WELL 

PRODUCTION FORMATIONl: A-JENKINS PRODUCTION FORMATION 2: 

IPCOF® 6939.000 MGF IPCOF® 

IPF® 0.000 MCF IPF® 0.000 MCF 

INTERVAL!® 3177-3187 INTERVALl: 0- 0 

INTERVAL2® 0-0 INTERVAL2: 0- 0 

ACID 0 GALLONS ACID® 0 GALLONS | 

SAND= 0 # SAND® 0 If 

GALLONS® 0 GALLONS® 0 

FRAC PRODUCTION® 0 MCFGPD FRAC PRODUCTION® 0 MCFGPd] 

• '.U 

HO. 1 CUMULATIVE PRODUCTION+ 6939.000 CODE®;^] 


Figure 15. Format for rectangular window computer 
output. 
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Analysis of the well stimulation data provides little 
evidence to support the hypothesis that surface linears can 
be related to subsurface reservoir fracturing. Although 
approximately 40 percent of the producing wells in the 
entire study area are hydraulically fractured, the average 
open flow potential for these wells is not significantly 
higher than the average for non-stimulated wells. Within 
the four separate areas selected for detailed analysis 
(Table 12), there appears to be no relationship between the 
open flow potential of hydraulically fractured wells and 
linear density, nor can a relationship be established between 
the open flow potential of stimulated wells and linear 
proximity. ^ 

B. Northern Region 

In parts of the northern region well control and linear 
density contrast are sufficient for use in examining the 
relationship between well yield and linear trends. To 
facilitate computer analysis, the northern one third of the 
test site was subdivided into half ‘-quad (7,5' quad maps) 
areas (Fig. 16). In addition, the various gas-bearing zones 
within the stratigraphic column were consolidated into 
manageable groups based on their stratigraphic position and 
reservoir and lithologic similarities: Upper Atoka, Lower 

Atoka, Sandstone Group, and Carbonate Group. 


Tafele 12,. Comparison of linear areas 



9415D0 56.2 sq/ km 












PRODUCTION PLOT 
LEGEND 


+- NON-PRODUCING WELL 
A -< 2 MILLION CUBIC FEET 

0=2 <5 MILLION CUBIC FEET 
O = 5 <15 MILLION CUBIC FEET 
□ = 15 MILLION CUBIC FEET OR MORE 

NUMBER OF PRODUCING WELLS = 1059 


Figure 16* Half-quad subdivision, northern region. 

Linears from Landsat winter scene band 7 
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The Atoka Forinationr which accounts for a large part 
of the gas production r contains a laterally persistent shale 
horizon between the Vernon and Sells sands. This shale is 
chosen as a boundary between the Upper Atoka and the Lower 
Atoka (LATOKA) , The upper Atoka units are very thin or 
absent in the northern margin of the test site and therefore 
are omitted from the analysis. The sandstone units in the 
Lower Atoka are: 


Lower Atoka (LATOKA) 

Sells Cecil Spiro 

Jenkins Patterson 

Dunn C Spiro 

Paul Barton 


The remaining productive horizons can be grouped by 
lithology and demarcated by unconformities. In descending 
stratigraphic order, the next group is predominantly 
sandstone and is referred to as the Sandstone Group (SAND ST). 
The Sandstone Group, bounded by unconformities, is from 
Pennsylvanian to Mississippian in age. It includes the 
following units 

Sandstone Group (SAND ST) 

Kessler Upper Hale 

Cannon Middle Hale 

Pitkin Lower Hale 

Wedington 

The Carbonate Group (GARB), the lowermost group, is 
predominantly limestone and dolomite and ranges in age from 


Mississippian to Ordoviciai). Tho Carbonate Group is also 
bounded by unconformities and includes the following units. 

Carbonate Group (CARB) 

Boone Everton-Sirapson 

Penters Arbuckle 

Hunton 

Another parameter, referred to as Total (TOTAL), is needed 
to represent the total production from the three groups for 
any certain well. 

1. Density Analysis 

For the density values, each linear within a half-quad 
was measured and the lengths were added to obtain a total 
linear length for the half-quad. This value is referred to 
as DENSl. Linears (or any part thereof) within a half-quad 
are hypothesized to have an effect on gas production and 
linears outside the half -quad are not. The total number 
of linears within the half -quad was calculated and is 
referred to as DENS2 (Table 13). 

By means of the Statistical Analysis System (SAS) (Barr 
et al., 1976), the Pearson correlation coeffieicnt test is 
used for a relationship between the production values and 
the density functions. 

DENSl - TOTAL 
DENSl - LATOKA 
DENSl - SANDST 
DENSl - CARS 
DENS2 - TOTAL ^ ^ ^ 


Table 13. 

llalf-guads and density measurements 

• 

11a If -Quad 

DEN SI 

DENS2 

Number 

(Total Length-Meters) 

(Total Number) 

1 

35305 

13 

2 

26890 

17 

3 

. 42409 

29 

4 

48445 

35 

5 

45854 

25 

6 

29268 

17 

7 

39909 

21 

8 

43689 

26 

9 

52370 

38 

10 

44055 

30 

11 

40274 

25 

12 

44024 

31 

13 

38994 

22 

14 

32896 

. 18 

15 

19878 

13 

16 

42226 

23 

17 

43567 

20 

18 

60396 

48 

19 

36189 

19 

20 

32805 

22 

21 

26402 

14 

22 

58079 

47 

23 

44726 

38 

24 

50732 

29 

25 

51707 

34 

26 

44695 

30 

27 

49695 

35 

28 

38415 

22 

29 

16311 

6 

30 

19421 

9 

31 

18201 

8 

32 

29085 

18 

33 

7683 

7 

34 

17591 

9 

35 

27530 

18 

36 

23994 

16 

37 

31707 

21 

38 

41524 

27 

39 

45091 

36 

40 

37896 

28 

41 

43537 

32 

42 

38963 

25 
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DENS2 - LATOKA 
DENS2 - SANDST 
DEWS2 - CARD 

If a half-quad has less than two wells it is not 
considered in the correlation analysis. Fourteen half-quads 
have two or more wells (table 14). 

A correlation coefficient (p) and significance 
probability level (a) were generated for each combination 
of the six parameters (Table 15), The null hypothesis states: 
•'There is no correlation between linear density and well 
production. " A correlation is assumed to be present if a 
relationship has a significance probability level greater 
than 95 percent (a < .05) or a correlation coefficient greater 
than 0.5 (p > 0.5). 

From Table 15, one significant correlationis found. 
Production values from the Carbonate Group (GARB) show a 
positive correlation (ot 0.0031, p = 0.72866) with the 
number of linear s per half-quad (DEMS2). Most Carbonate Group 
production is from the relatively low porosity, impermeable 
Boone Formation (Van Den Heuval, 1979) and therefore may 
depend on the number of fractures. 

* - 

, A relationship is suggested (a =0.0756, p “0.51827) 
between the Sandstone Group (SANDST) and the number of linears 
per half-quad (DENS2) . Dnlike the units in the Lower Atoka, 
the Sandstone Group includes more mature, calcarGOus 
sandstone. The calcareous cement of this sandstone could 
increase its responsiveness to natural fracturing. 
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coefficient ^ density. 
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There is also a positive correlation (a = 0.0253, 
p » 0.59359) between total production (TOTAL) and total 
nurober of liner as per half -quad (DENS2) . This finding is 
of questionable value because the Garhonate Group production 
ma)ccs up a large part of the total production. No other 
relationship of significance are found. 

2. Proximity Analysis 

For this analysis, two distance parameters are used. 

One is the di^^tahce from a given well to the nearest linear 
(DISTl) and the other is the distance from a given well to 
the nearest linear intersection (DIST2) (Table 16) . Two 
linears cire not considered intersecting if they are more 
than two line widths (0.5 mm) apart. 

The definitions for both distances are based on the 
concept of an expanding circle with a given well as the 
center. As the circle expands it eventually >meets a linear 
and a linear intersection. The distance is the radius of 
the circle at the point of contact with the linear and linear 
intersection. 

By means of the Statistical Analysis System (SAS) 

(Barr et al. , 1976) the Pearson correlation coefficient is 
used to test for a relationship between; 

DISTl - TOTAL 
DISTl - LATOKA 
Distl - SANDST 
DISTl - GARB 
DIST2 - TOTAL 


Table 16. Data used for Pearson correlation coefficient 
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Table 16. Continued. 




DIST2 - LATOKA 
DIST2 - SANDST 
DIST2 - CARB 

As before, a correlation coefficient (p) and significance 
probability level (a) are generated for each combination of 
the six parameters (Table 17). The null hypothesis states: 
"There is no correlation between distance to linear (and 
linear intersection) and well production," A correlation 
is assumed to be present if a relationship has a significance 
probability level greater than 95 percent (a < .05) or a 
correlation coefficient greater than 0.5 (p>0.5). 

From Table 17 a negative correlation (a = 0.0450, 
p = -0.253) is suggested between the Carbonate Group (CARB) 
productionand the distance ^o the nearest linear intersection 
(D1ST2). In other words, the closer a well producing from 
the Carbonate Group is to a linear intersection, the greater 
its production rate. Most of the producing wells in the / 
Carbonate Group are 305 to 2134 m from the nearest linear 
intersection and wells beyond 2134 m are dry holes. The fact 
that many dry holes are also between 305 and 2134 m from the 
nearest linear intersection explains the marginal correlation 
(a > 0.0450, p= -0.25348) . 

The important relationship is between production from 
the Carbonate Group (CARB) and the numl>®J^ of linears per 
half-quad (DENS2 ) . The marginal correlation of Table 17 
reflects this relationship. As the number of linears per 
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Table 17. Pearson correlation coefficient - distance 
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half-quad increases, the frequency of intersecting linears 
also increases. In other words, if a producing well is in 
an area of high linear densi-ty, there is probably a linear 
intei^section within 2134 ra. This relationship has application 
for locating both exploration and development wells, in 
sandstone and carbonate strata, the distance to the nearest 
linear or linear intersection can be disregarded. A well 
is wore likely to be productive if it is in an area of high 
linear density ai*.d if it penetrates structurally elevated 
carbonate strata. 

C. Practical Application 

Within the northern region, a relatively large number of 
producing gas wells are clustered in the vicinity of half- 
quads 34 and 22 (Fig. 16) . Visual analysis of Figure 16 
indicates a contrasting relationship for linears in these 
same quads . Whereas half-quad 34 (immediate vicinity of 
Rock Creek gas field) is characterized by few linears, 
half-quad 22 has a relatively large- number of linears, 
especially in a north-south orientation. In an attempt to 
define the relationship between linears, subsurface structure, 
and gas production in these two contrasting areas, subsurface 
structure maps were constructed from mechanical well logs 
(induction electric, formation density, and gamma ray/neutron ) 

1. Batson Gas Field 

Although a number of wells produce gas from sandstone 
units within the Atoka Formation, the Batson Field is the 
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only area of significant Boone Formation production in the 
entire Arkoina basin. Figure 17 is an enlarged plot of 
linears and gas-producing wells in all rock units in the 
Batson field area. 

Wells producing gas^from the Boone Pormation are on 
structural highs in the Batson field , and wells not on 
structural highs have produced water. From Figure 18, the 
-900 foot (below sea level) contour was chosen to distinguish 
wells that are structurally high from those that are 
structurally low. There is a positive correlation between 
production from the Boone Pormation and well location with 
respect to structure, with two exceptions . The well in SW 
1/4 of Sec. 32 , T. 12 N. , R. 24 W. (Pig. 18) is a dry hole 
in the Boone Formation in spite of its structurally high 
position. Examination of subsurface cross-sections 
constructed from mechanical well logs indicates that the ; 
well intercepts a fault, which appears to eliminate any 
production possibility. The well in Sec. 29, T. 11 N. , R. 

24 W. (Pig. 18) is a producing well but is in a structurally 
low area of the Boone Pormation. The well is outside the 
area where the contours can be drawn with confidence (limited 
control points). Moreover , production records show this 
particular well to be a very low producer (72.4 MCFGPD) . 

2. Rock Creek Field 

Figure 19 is an enlarged plot of linears and gas-producing 
wells in the Rock Creek field. A structure contour map was 


81 













Non - p ro d u c i n g we I I 


o 


H O 
H (U 

w ? U 

« 

Z 

c w o 

£ w « 

c >t 


W ^ 

m o 

ig 0i 

e a 

O *H 
M 4J 

*W 'r( 

c 

fd >H 

o 

■H 

o 

(U -H 


flj 4J »0 
0) d H 
COO) 
•H 0 -H 
>4 Hm 


c O 


V 

V CVl 


+ < ♦ 


conatructed for the top of the Boone Formation (Pig. 20). 
Several structural highs are present but# unlike veils in 
the Batson field# wells on these highs do not produce from 
the Boone Formation. In the entire area only two wells 
produce from the Boone Formation (Pig. 20) and both are 
relatively low producers. Inspection of the samples from 
the productive horizon of the well in Sec. 24# T. 11 N# K, 

20 W. (Pig. 20) shows minor fracture porosity. For the 
producing well in Sec. 22# T. 11 N.# R. 27 W. (Fig. 20)# 
fracture porosity is not found but production records show 
the well to be a very low producer (80 MCPGpd) . 

3. Data Analysis • 

The Boone Formation (up to 60-70 percent chert) yields 
gas in the Batson field which is situated on a structural 
high coincident with high linear density. However# in the 
Rock Creek field# the Boone Formation is on a structural 
high in an area of low linear density uud is essentially 
nonproductive . A study of mechanic; al logs from both fields 
indicates that the Boone interval is generally lithologically 
consistent throughout. 

The Pearson correlation coefficient indicates a positive 
correlation between the number of liners in a given area 
(linear density) and production from the Boone Formation; 
however , there is no significant correlation between production 
and the distance to the nearest linear. Because the Boone 
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Structure contour map for the top of the Boone 
Formation, Rock Creek field . 
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Formation is a brittle# high-chert-content limestone with 
very low porosity and permeability (Van Den Heuval# 1979)# 
it is reasonable that the high surface linear density is a 
reflection of fracture porosity/permeability in the subsurface 
Previous studies in other areas (Harris et al.# 1960; 
Stearns and Friedman# 1972; Gorham et al.# 1979) indicate 
that the strongest surface fracturing is on the flanks of 
folds where curvature is at a maximum. In addition# Gorham 
et al. (1979) sr.ggest that open fractures had to develop 
where tensional joints form at places of maximum curvature 
of beds (i.e.# where there is greatest rate of change of 
'dip# not necessarily where the dips are steepest). The 
Batson field is on the hinge line of a regional east-west- 
trending monocline resulting from the transition from 
essentially horizontal strata of the northern Arkansas 
structural platform to the southward-dipping strata of the • 
Arkoma basin, open fractures are believed to develop where 
tensional joints form at places of maximum curvature (hinge 
line) of the beds (Gorham et al., 1979). Figure 17 does not 
indicate a prominent east-v;est-trending linear pattern. 
However# if the fracture-development hypothesis is correct 
and such an east-West-trending pattern is present in the 
subsurface rocks# the prominent north to northeasterly linear 
pattern defined from Figure 17 would be superimposed upon it. 
The somewhat orthogonal nature Cf such a pattern appears to 
be supported by field work# as a general north-south and 


east“*wGst pattojrn of surfape joint data has boon recorded 
(Figure 4 ) . 

As discussed in Section XVf the cast-west-trcnding non** 
systomatic surface joints are probably not charactoristically 
indicative of tensional stress conditions; however, because 
of the structural position of the Batson field, relatively 
tight fractures way have been changed to open fractures. 

I'he structural position of the Batson field on a regional 
n\onocline, coincident with a local structural high with 
inareased linear density, provides empirical support for an 
intorrelationship between remote sensor-derived linears 
and fracture production. The Rock Greek field, although on 
a local structural high, is off the hinge line of the 
monocline in an area of low linear density, and has minimal 
evidence of fracture production. These findings have direct 
application to exploration and development of the Boone 
Pormation in the northern Arkoma basin. More important, 
the findings agree with those of Gorham et al. (1979) that 
exploration efforts for traps with fractured reservoirs should 
focus on areas of maximum curvature of folds, especially 
monoclines. 

VII. GEOPHYSICAL STUDIES 
A. Gravity and Magnetic Surveys - 

Gravity and magnetic surveys were completed along the 
southern extension of three prominent northeast- trending 
lineaments (Drakes Creek, Bonca and Lurton) which appear to 


transect several structural domains (Fig. 21). These 
lineaments appear to extend from the northern Arkansas struc'- 
tural platform southward into the Arkansas part of the Arkoma 
basin. North of the Arkoma bas*n parts of these lineaments 

have been mapped as normal faults that are downthrown on the 

* 

southeast. Although there is no identifiable surface evidence 
of faulting for these lineaments in the Arkoma basin, geophysi>- 
cal studies were conducted in an attempt to establish their 
southward extension into the gas-bearing strata of the Arkoma 
basin. 

1. Location and Data Collection 

Survey patterns were established in two separatG areas 
(Areas I and tt, Fig. 21) for collecting gravity and magnetic 
data. Gravity data were obtained with a Worden gravimeter 
and standard error corrections were applied before values 
were plotted on base maps. Magnetic data were obtained with 
a Varian proton procession magnetometer. Short-term variations 
in the earth's magnetic field were established approximate ly 
every two hours to increase the accuracy of readings at 
individual stations.. In addition to contouring data, gravity 
and magnetic maps were evaluated with the aid of a Fortran 
IV-llenderson-2 program (Henderson, 1960) . This program 
generates both first and second derivative maps and upward 
and downward continuation maps. 








2. Data Analysis 

In the northern part of Area I, a small gravity trough 
with a maKimum magnitude of 3 mgals was delineated more or 
less on trend with the Drakes Creek lineament. Associated 
with this feature is a fairly uniform decrease in magnetic 
intensity/ sloping to the southeast. It is possibly 
indicative of displacement of basement downward on the 
southeast along the Drakes Creek lineament. In the southern 


central part of Area 1 is a roughly circular magnetic anomaly 
of at least 100 gammas magnitude (circle in Area I , Fig. 21). 
This anomaly is adjacent to the possible southwest extension 
of the Drakes Creek lineament. Modeling led to the 
conclusion that the source of this anomaly is either an 
erosional remnant or possibly a rhyolitic extrusive pile on 
basement (Watkins / 1978; Schumacher, 1979) . 

Magnetic profile data across the sou thv;e stern extension 

f 

of the Ponca lineament in Area II are somewhat similar to 


profile data obtained by Smith (1977 ) along the Ponca 
lineament's northeastemmost extension (Area III, Fig. 21) . 
Smith (,1977) noted that the Ponca lineament (Area III) 
coincides with the linear arrangement of lead-zinc mines in 
the Ponca-Boxley district. On the basis of gravity and 
magnetic data analysis. Smith suggested that the northern 
part of the Ponca lineament is probably a first-order fracture 
and/or shear zone that penetrates both Paleozoic and 
Precambrian rock. By inference, then, the Ponca lineament may 


pirovide a continuous shear zone from Area III through 
Area II. 

A large, nearly circular magnetic high is centered 
about 1,6 km north of Clarksville, Arkansas, intersecting 
the trend of the burton lineament (circ].o in Area II, Pig, 

21) . This anomaly covers nearly three townships and has a 
magnitude of approMimately 450 gammas above the regional 
field (McBride, 1980), Preliminary modeling suggests that 
an intrusivo body with relatively high magnotic susceptibility 
is present in the Precambrian basement. Because of the 
masking effect of this large anomaly, the subsurface extent 
or presence of the southwestern extension of the burton 
lineament cannot be determined. 

B. Results 

Gravity and magnetic data analysis provides evidence to 
support the southwestward continuation of the Drakes Creek 
lineament into the gas -bearing strata of the Arkoma basin 
(Area I) . The extension of the Ponca lineament (Area II) 
may coincido with inferred magnetic gradients related to 
basement shear zones, but similar evidence is lacking to 
support the presence of the burton lineament. It sliould be 
noted, however , that parallelism between northeast-southwest- 
trending lineaments and magnetic gradient anomalies has 
been reported in the next physiographic province on the east , 
the northern Mississippi embayment (O'beary and Simpson, 
1977) . In addition, recent analysis of gravity and 
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aeromagnotic data ftrom the notthorn part of the Mississippi 
entbayment suggests that basement tectonics are largeXy 
responsible for the close association between geophysical 
anomalies and intrusions of basic and ultrabasic rocks 
(Braile et al.# 1979; HendrickSi 1979; Hildenbrand et al,, 

1979) . 

The relationship between lineament extension, gas 
production, and reservoir characteristics is not apparent. 

There is no correlation between gas production and proximity 
to the three lineaments extended into the Arkema basin. 

VIII. SEASAT IMAGERY 

Seasat swath Rev 795 provided a 100-km fwath of SAR 
imagery across southern and western Arkansas. Unfortunately, 
the digitally processed SAR covers only the northwesternmost 
corner of the test site area. No comparison with linear 
data sets could be accomplished. 

IX. SUMMARY AND CONCLUSIONS 

Surface joints in the Arkoma basin consist of two 
orhtogonally oriented sets. One set (systematic) , trending 
N10®W, consists of straight, planar, through-going j oints . 
Plumose structures on many of these joints suggest a ten s ion a 1 
origin . The other set of joints (non-systematic) has an 
average orientation of N86®E and consists of short and often 
curved joints which cut the rock between the systematic joints. 
Along the northern flank of the Arkoma basin, where gas 



* 


pjTocluction appears to be from a fractured reservoir, the 
north-south-trending joint set is somewhat coincident with 
a prominent north-south linear trend, in general, the 
influence of surface jointing on linear s throughout the 
Arkoma basin is not obvious. 

Of Landsat bands, band 7 detects the greatest number of 
topographic linears. However, the variation among the 
three operators is least with band 6. The wintor liandsat 
scene consistently provides for the detection of more linears 
than the summer handsat scene because of the lower solar 
elevation during the winter. The uniform distribution stretch 
enables interpreters to detect more linears than the Gaussian 
stretch. The Gaussian stretch enhancement decreases contrast 
in the middle-gray range, thereby suppressing linear 
expression. In summary, the best Landsat combination is the 
winter scene, band 7, uniform distribution stretch. 

Of the SAR data products, the VH SAR radar mosaic 
provided for detection of the most linears; however , none 
or the SAR enhancement products is significantly better than 
the others. Two of three operators detected fewer linears 
from the radar/Landsat merge than from the Landsat data set 
alone; however, because of operator variability, the results 
are inconclusive. A black and white band 7 wintor Landsat 
scene merged with radar mosaic (VH) , having a look-direction 
orthogonal to the Landsat solar azimuth, is suggested as a 
complementary image merge . Both of these data sets alone 
allowed for the detection of most linears. 
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Radar look-direction and Iiandsat solar azimuth illumin- 
ation direction tend to suppreRS topographic linears more 
parallel with illumination direction and to accentuate those 
having a more perpendicular relationship. 

The operator variability studies determined that 

* 

generally the frequency, average length, and coincidence of 
linears differ greatly among operators, although major linear 
orientation trends show some degree of similarity. This 
observation has also been made by Podwysocki (1975) and 
Siegal (1977). in general, linear analyses tend to be 
subjective and highly variable. Even when the operator 
variables are minimized, variations between remote sensor 
data sets may produce significant differences in interpre- 
tations. No single technique of interpretation or single 
type of imagery seems to be universally best for minimizing 
operator variability. However , maps derived from linear 
density analysis tend to be somewhat less influenced by 
operator variability. Geneo.ally, conservative linear "mappers" 
tend to infer fewer linears per unit area, but relative 
concentrations are sometimes agreed upon among operators . 

Gravity and magnetic data analysis provides evidence to 
support the southward continuation of two of three major 
lineaments into the gas-bearing strata of the Arkoma basin. 

The relationship between lineament extension, gas production, 
and reservoir characteristics is not apparent. No correlation 
is found between gas production and proximity to the three 
lineaments extended into the Arkoma basin. 
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The southern region oi? the Ark<?^sas Arkoma basin proper 
is characterized by high well density but few linears. 

Little relationship is discernible between surface structure 
and gas production# and no correlation is found betwetin gas 
productivity and linear proKimity or linear density as 
deterinined from remote sensor data. Stratigraphic gas 
entrapment is most probably related to depositional environ- 
ment# and is the main factor influencing porosity and 
permeability in the predominantly sandstone reservoirs. 

in the northern Arkoma basin a positive correlation 
is found between the number of linears in a given area and 
production from clierty carbonate strata. Most such production 
is from the Boone Formation. The high chert content of the 
Boone Formation makes it particularly susceptible to fracturing 
and the development of fracture permeability. No correlation 
is evident between linear density and any of the other 
gas~bearing units. These other gas-bearing units are 
predominantly sandstone# and are somewhat similar lithologically 
to the sandstone reservoirs in the southern Arkoma basin. 

Linear dettsity analysis should not be the single basis 
for petroleum exploration on the northern flank of the Arkoma 
basin. Regional and local structural position and surface 
joint orientation are equally important data sets to be 
evaluated in an exploration program. However # the probability 
of economic savings in the exploration for fractured reservoirs 
(similar to the Boone Formation) in other parts of the world 
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can be considorablo if lineaj? density maps derived from 
radar/Landsat merges » are Used to target prospects for 
further goologicai and geophysical inv^sistlgations, 

X. RECOMMENDATIONS 

The following recoirmi'endations are based on the experience 
gained in the performance of this contract, 

1, Simultaneous multifreguency-multipolarization radar 
imagery should be obtained in an area where linear s are 
numerous, but whore operator variability can be held to a 
minimum. The test site should be flov/n in such a fashion 
to provide four orthogonal looX-directions, and at least 
three contrasting depression angles (20®, 45®, 70®). This 
data set would proviiv*' nufficient imagery to allov/ quantitative 
analysis. I 

■ : I ' ■ is 

2 . Individual interpretations of identical imagery are 
highly variable and tend to be subjective. The applicability 
of using machine |brocessing needs to be determined, 

3. The optiritum landsat data set (winter scene band-7) 
and L-bandSAR imagery should be merged and evaluated with 
respect to the radar/Landsat combinations reported on in 
this study. 

4. Space Shuttle radar imagery (with a look-direction 
different from Landsat solar azimuth preferably^ orthogonal) 
should be viewed in a complementary role for topographic 
linear analysis where Landsat imagery is available. 


5. Linear density rather than linear proximity should 
be emphasized in the search for fractured reservoirs; 
however f linear analysis should be used in concert with 
Other geologli? and geophysical information. 
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appendix 


Index to Linear Histograms 

Histogram 

Number 


Landsat Winter, Gaussian Stretch 

1 . Band 4 

2. Band 5 

3. Betnd 6 

4. Band 7 

Landsat Winter, Uniform Distribution Stretch 

5. Band 4 

6. Band 5 

7 . Band 6 

8. Band 7 

Landsat Summer, Gaussian Stretch 

9. Band 4 • 

10. Band 5 

11. Band G 

12. Band 7 

Landsat Summer, Uniform Distribution Stretch 

13. Band 4 

14 . Band 5 

15. Band 6 

16. Band 7 

17 . Landsat Winter, Enhanced Color 

18. Landsat Winter, Standard Color 

19. Landsat Summer , Enhanced Color 

20. Landsat Summer, Standard Color 

SAR, L-Band Imagery — May 1977 

21. Uncontrolled Mosaic (HH) 

22 . Uncontrolled Mosaic (HV) 

SAR, L-Band Imagery—April, 1978 HH Mosaic 

23. Response Range Dependence Removed^^^^^^ 

24. Uniform Distribution 

25. High Pass Filter 

26. High Pass, 20% Added Back^^^^^^^ ^ 
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SAU, L-Dand Imagery-'~ApriX| 1978 VII Mosaic 

27. Response Range Dependence Removed 

SAR Radar/Iiandsat Composite 

28. SARHin and Dandsat Winter, Enhanced Color 

29. SAR'^HH and Eandsat Svimmer, Enhanced Color 

Modified Area Histogram Coverage 

30. handset Summer, Enhanced Color 

31. handsat Winter, Enhanced Color 

32. HH Mosaic 
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